The role of the Berry Phase in Dynamical Jahn- Teller Systems 



Nicola Manini 
European Synchrotron Radiation Facility, B.P. 220, 



Grenoble Cedex, France 



oo 



Paolo De Los Rios 

Institut de Physique Theorique, Universite de Fribourg, 1700-CH Fribourg, Switzerland 

(February 1, 2008) 

The presence/absence of a Berry phase depends on the topology of the manifold of dynamical Jahn- 
Teller potential minima. We describe in detail the relation between these topological properties 
and the way the lowest two adiabatic potential surfaces get locally degenerate. We illustrate our 
arguments through spherical generalizations of the linear T ® h and H ® h cases, relevant for the 
physics of fuUerene ions. Our analysis allows us to classify all the spherical Jahn- Teller systems 
with respect to the Berry phase. Its absence can, but does not necessarily, lead to a nondegenerate 
ground state. 
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I. INTRODUCTION 

The traditional field of degenerate electron-lattice in- 
teractions (Jahn- Teller effect) in molecules and impurity 
centers in solids has attracted new interest in recent 
years, excited by the realization of new systems which call 
for a revision of a number of commonly accepted beliefs. 
A whole range of icosahedral molecular systems including 
Ceo ions and some higher fuUerenes, thanks to the rich 
structure of the symmetry group, are characterized by 
up to fivefold-degenerate representations of the electronic 
and vibrational states of the isolated molecule/ion. New 
Jahn- Teller (JT) systems have therefore been considered 
theoretically, |||-^ with intriguing features. Q-Q] A par- 
ticularly surprising property has been demonstrated re- 
cently: the possibility of a symmetry switch of the molec- 
ular ground state [^|J|] , connected to the absence of Berry 
phase in the coupled dynamics of electrons and vibra- 
tions. 

As it is well known, the molecular symmetry, reduced 
in the static JT effect with the splitting of the electronic- 
state degeneracy, is restored when the coherent tunneling 
between equivalent distortions is considered, in the dy- 
namical Jahn- Teller (DJT) effect. In this context it was 
commonly accepted an empirical "symmetry conserva- 
tion rule", sometimes referred to as "Ham's theorem", 
stating that the symmetry of the vibronic DJT ground 
state, at all coupling strengths, remains the same as that 
of the electronic multiplet prior to coupling: ||^ all lin- 
ear JT systems known till a few years ago, for single- 
electron occupancy, systematically satisfy this empiric 
rule. It was understood recently that this phenomenon, 
not automatically implied by the DJT physics, is in re- 
ality a fingerprint of a Berry phas e [p| in the entangled 
electronic-phononic dynamics. P,|7|,p^|~Consequently, this 
geometrical phase seemed a universal feature of the DJT 
systems. 



It came unexpected the discovery of the first dynamical 
JT system without Berry phase, showing a nondegener- 
ate ground state in the strong-coupling limit. ||^,^ This 
is the model that in spherical symmetry is indicated as 
2?*^^) ® d^'^\ where electrons of angular momentum L = 2 
interact with vibrations also belonging to an / = 2 rep- 
resentation. This system is relevant to the physics of 
fuUerene ions Cqq, where the 5-fold degenerate electronic 
state has iJ„ icosahedral label and the quadrupolar dis- 
tortions correspond to some of the kg modes. It has 
been shown both analytically and numerically that, for 
increasing coupling, a nondegenerate state in the vibronic 
spectrum moves down, to cross the 5-fold ground state 
at some finite value of the coupling parameter, thus be- 
coming the ground state at strong couphng 0,^. 

In this work, we review the mechanism of the Berry 
phase in degenerate electron- vibration coupled systems, 
and uncover the detailed reason of the absence of this 
geometric phase in the P^^-' (8) rf^^-* system, hence drawing 
some natural generalizations. We structure the paper as 
follows. In Sect. ||, we introduce the basic JT Hamil- 
tonian, along with the fundamental transformations to 
recast it in a semi-classical form. Section [II reviews 
the mechanism responsible for the presence of the Berry 
phase in most DJT systems, and its consequences for 
the low-energy part of their spectrum. In Sect. IV , the 
mechanism allowing to get rid of the Berry phase in some 
systems is unveiled, and a class of models sharing this 
property is proposed. Section ^ reports some numerical 
results confirming and complementing the predictions of 
Sect. In Sect. VI, we discuss possible generalizations 
of our model to higher-order couplings and discrete sym- 
metry groups. Finally, conclusions are drawn in Sect. 



VII 
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II. THE MODELS 

According to the general theory of the JT effect, ||] 
an TV-fold degenerate electronic level corresponding to a 
representation T of the molecular symmetry point group 
interacts with all the vibrational modes corresponding 
to representations {A} contained in the symmetric part 
of the direct product of F with itself. For simplicity, in 
this work we restrict to only one degenerate mode. The 
Hamiltonian for this "F (g) A" case reads: 



|A| 



(1) 



where &J / bi are the creation/destruction operators for 
the harmonic vibrational mode component i, and -ffo-v 
is the interaction Hamiltonian, which, to linear order in 
the boson operators, writes 



|A| |r| 



J2 [hc]ck{Ai\Tjrk) + h.c. 



i—1 j.k—1 



(2) 



Cj is the fermion operator for orbital j, {Ai\TjTk) are 
the Clebsch-Gordan coefficients of the symmetry group 
of the problem, and g is the dimensionless coupling pa- 
rameter. For the purpose of illustrating our analysis, it 
is particularly convenient to stick to the case in which 
the symmetry group is that of three-dimensional rota- 
tions, S0(3). In the following, therefore, we focus on the 
coupling of an electronic state of angular momentum L 
(whose representation F is indicated as P^^^) to a degen- 
erate vibration of angular momentum I (of representation 
A = For this case, the Hamiltonian reads 



m— — l 



k.k' = -L 



[6L + (-ir&-r„]4cfc,(/m|ifc,L-fc') } , (3) 

where some of the symmetries of the 5*0(3) Clebsch- 
Gordan coefficients are implied. 

The Hamiltonian (^) is suitable for perturbation calcu- 
lations (small g values) and as a starting point for numer- 
ical diagonalization methods, such as the Lanczos tech- 
nique. Yet, the Berry phase is a semiclassical concept, 
useful only in the medium/large g regimes. Here it is 
convenient to switch to a real representation of the vi- 
brational degrees of freedom. 

We apply two unitary transformations, both on the 
electronic operators, and the vibrational ones. We define 
a new set of electronic operators, Cm (and consequently 
their hermitian conjugate c}„), m = —L,...,L via the 
transformation 



Co = Co 
\ cxp 



(4) 



V2 



I —I 
1 1 



, m > 



The second transformation expresses the 2Z -I- 1 boson op- 
erators Bm = ^Jri + ("1)™ ^-m in tcrms of the Hermitcan 
"coordinate" operators as follows: 



Bo = V2go 



(5) 



i 1 



q- 



771 > 



The remaining components foj^^ — (— l)™6_m can be ex- 
pressed correspondingly in terms of the momentum oper- 
ators pm conjugate to qm. We stress that these transfor- 
mations are by no means unique. Any further orthogonal 
transformation of the real coordinates, leads to equiva- 
lent results. 

Eventually, the Hamiltonian operator (|^) transforms 
into the form 



1. 



I 



(6) 



with 



1 



m—-~l 



L 

^ E 



(7) 

M) 
For 

brevity, we introduce the symbols 7^^. = (/ j+k\L j, L k), 



It is straightforward to compute the (2L -I- 1) x {2L 
coupling matrices V^"^\ for any value of L and 1. 



are: 



0. In this notation, the matrices for the L = 2 case 
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Note that the forms of V'^^\ and ^(±2) ^^pp^y 

to the coupUng to both I = 2 and / = 4 vibrons (of 
course, the numerical values of the coefficients are dif- 
ferent), while the V'^^'^^ and V^^^^ matrices are relevant 
only for 1 = 4. 

For larger L's, the structure of the central 5x5 block 
(corresponding to indexes j = —2 to 2) of the corre- 
sponding coupling matrices is conserved, additional ma- 
trix elements being added externally in a simple way. For 
example, for any L and I, the V^^"^ matrix is diagonal with 
matrix elements F^f'^'' = (— 1)'^7^!'_^.- Similarly, \il = 2-L 
then the only nonzero elements of the matrices are 



V_ 



-i^andKi-'^ 



±i,±i ~ ^ "-"-^"^ " ±1,^1 ~' j^- 
As anticipated, the form (^,0) of the Hamiltonian is 
suitable to a semiclassical treatment, in the spirit of the 
Born-Oppenheimer (BO) approximation, i.e. a factoriza- 
tion of the electronic "fast" dynamics from the "slow" 
distortions which are quantized as a second step. In 
this scheme, we briefly review the Berry phase mecha- 
nism, and its consequences on the factorized dynamics. 



III. THE BERRY PHASE 

The traditional BO scheme assumes that the elec- 
trons, moving much faster than the ions, follow adia- 
batically the ionic (vibrational) motion with the only 
effect of generating a potential energy for the vibra- 
tional motion. This approximation relies on separa- 
tions between consecutive electronic levels that are much 
larger than the typical vibrational energies fiio. In a 
JT problem, the BO treatment starts with the diago- 
nalization of (^ in the electronic degenerate space, at 
fixed distortion field i.e. the diagonalization of the ma- 
trix S = ^qmV'^™'^ ■ Each electronic eigenvector ji/;^) 
of S, of eigenvalue A^, generates a BO potential sheet 

14(g) = huj ^Em^m + ' including the har- 

monic potential of the free vibrations. At strong enough 
coupling gr, the separation of the potential sheets becomes 
large enough that the adiabatic motion can be assumed 
to always follow the lowest BO potential sheet, while vir- 
tual electronic excitations may be treated as a small cor- 
rection. For Hamiltonians of type (||), the set of points 
9min of minimum potential energy, i.e. the classical stable 
configurations, constitutes a continuous manifold, often 
denominated Jahn- Teller manifold (JTM), and the value 
of the lowest BO potential there is the classical JT sta- 
bihzation energy Scias = Kiin(gmin)- 



On the other side, due to time-reversal invariance of , 
the space of all possible (normalized) electronic eigen- 
states can be represented by an (hyper-) sphere in the 
[2L + l]-dimensional real space (see Fig. Ilj). The BO dy- 
namics realizes an adiabatic mapping of the vibrational 
space into this electronic space. Indeed, every point q 
on the JTM (in the vibrational space) is associated to the 
electronic wave function |'0min('7)), corresponding to the 
lowest eigenvalue An 
matrix S: 



of the electron- vibron interaction 



^{q) |'/'min(<z)) = A,nin(g)|V'min(g)) 



(13) 



This adiabatic mapping is two-valued, since opposite 
points ±|'0,„in(q)) on the electronic sphere give the same 
JT stabilization energy, thus corresponding to the same 
optimal distortion on the JTM. This identification of the 
antipodal points through the mapping is the mechanism 
allowing the JTM to have different connectedness jlj] 
from the electronic sphere. The latter is of course sim- 
ply connected, i.e. any closed path, or loop, on it can 
be smoothly contracted to a single point. The JTM, in- 
stead, may well be multiply connected, i.e. it can have 
intrinsic "holes" in its topology. In particular, in addi- 
tion to the contractable loops (such as that mapping on 
TTi on the electronic sphere sketched in Fig. |^), on the 
JTM we have the nontrivial class of those loops mapped 
on a path going from a point to its antipode on the elec- 
tronic sphere, such as 7r2 in Fig. |l|. In the traditional DJT 
systems {E®e, T®K), these two classes of paths are topo- 
logically distinct: loops belonging to class 2 may never 
be smoothly deformed into loops belonging to class 1, 
whence the multiple-connectedness property of the JTM. 
We see therefore that this multiple connectedness is in- 
timately related to the mapping between the JTM and 
the electronic sphere. The electronic sign change charac- 
terising the loops belonging to class 2 is a case of Berry 
phase [||. 

When also the vibrational motion is quantized, the 
overall (vibronic) BO wave function is factorized in the 
direct product of the electronic adiabatic state times the 
wave function for the slow degrees of freedom q: since the 
vibronic wave function is a regular, single-valued func- 
tion, the q degrees of freedom must cope with the elec- 
tronic phase change, which acts as a special boundary 
condition to quantization. As a consequence the motion 
on the JTM is constrained by special selection rules. For 
example the JTM of the simple E ® e system is a circle: 
the low-energy vibronic spectrum is indeed a spec- 
trum as for a circular rotor, but the Berry phase implies 



j = ±^,±|,..., instead of .7 = 0,±1,±2,... as for an or- 
dinary quantum rotor. |^,|l^ Similarly, the JTM of the 
T®h (i.e. 2?^^-' , in the spherical language) is equiv- 
alent to a sphere, [p|Jl^] but out of all the states, labeled 
by J, M, of a particle on a sphere, the Berry phase keeps 
only the odd-J ones. [^,^jlj] Note in particular that in 
these examples the presence of a Berry phase rules out 
a nondegenerate ground state. The same symmetry of 
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the degenerate electronic state, prior to the vibronic cou- 
pling/distortion, is enforced to the strong coupling DJT 
ground state. 



IV. ABSENCE OF THE BERRY PHASE 

As anticipated above, though a very frequent phe- 
nomenon, the Berry phase is not automatically implied 
by linear JT Hamiltonians (0). The above discussion 
should make clear that the absence of the Berry phase in 
a DJT system is linked to a mechanism leading to equiv- 
alence of the paths in the class 1 and 2. This mechanism 
should also coexist with the two-valued adiabatic map- 
ping sketched in the previous section. The solution of 
the riddle is provided by a point on the JTM where 
the mapping is degenenerate, i.e. it links qd not just to 
a pair of opposite points ilipminiQci)) on the electronic 
sphere, but to the whole circle of linear combinations 
cos 61 IV'i('fd)) +sin6' |V'2('fd)) of two degenerate electronic 
eigenstates (such as, for example, A in Fig. If such 
a point is present, any loop of class 2 on the JTM may 
be deformed smoothly so that its image on the electronic 
sphere becomes half this circle, thus the single point qd 
on the JTM: class 2 loops are therefore equivalent to class 
1 (contractable) loops, and, therefore, the JTM is simply 
connected. No Berry phase is possible in such a case. 

The possibility of such degenerate points may be con- 
sidered via careful analysis of the structure of the mul- 
tiple BO potential sheets given by the eigenvalues 
of the electron-vibration interaction operator. In most 
"classical" cases of Berry-phase-entangled linear JT sys- 
tems {E ^ e, T ^ h, ...), the lowest potential sheet re- 
mains separated from the next lowest one by a finite 
gap throughout the JTM. Conical intersections between 
these two sheets take place at some point in the distortion 
space, far from the potential minimum. In the E®e, for 
example, such a point is the origin q — Q. |jl^ There is, 
however, a second possibility: somewhere on the JTM, 
the lowest BO potential sheet gets tangentially degener- 
ate to the next lowest BO potential sheet. The energy 
difference between the two lowest levels is quadratic in 
the distance from the degenerate point. 

This possibility is indeed realized in our spherical 
model. Take for example the case L ^ 2, I ~ 2 and con- 
sider the point qd — (o, 0, —^9, 0, 0^ Simple inspec- 
tion of the coupling matrix S (||) shows that its lowest 

eigenvalue -\qd\l2-2 = -\ld\{l 0\L 2,L - 2) = is 
twofold degenerate. The classical BO potential energy is 
Ec\as = Vminiqd) = —jiQ^^'-^i and no other q yields lower 
potential energy than this: by definition, qd does belong 
to the JTM. By inspection of Eqs. (p[]lO|), it is possible 
to verify that the direct coupling between the two degen- 
erate electronic states (element 5±2,=f2 of the couphng 
matrix) vanishes identically: moving away from qd, these 
two states remain degenerate to first order, and the de- 
generacy is only lifted by indirect, second order coupling 



to the other states, qd is indeed a tangency point. With 
this, we have detailed the mechanism for the absence of 
a Berry phase in the I?*^^^ ® d^^^ model. 

It is natural to search possible extensions of this mech- 
anism to other cases of spherical DJT models (^. Eqs. 
(p|-|l^) show that on the qo axis, the electronic states are 
pairwise degenerate (except for the Cq|0) state). We are 
interested in the pair corresponding to the lowest and 
next lowest BO sheet, thus to the maximum numerical 
value of the Clebsch-Gordan |7;J; For the P'^) ^^(2) 
described above, the maximum is indeed ^2-2- 0] 
the contrary, for the V'^^^ (g) d^'*) case, the largest coefR- 

1 /2 

cient is 7,^0 = (M ) ' corresponding to the singlet state. 
Thus on the go axis there are no degeneracy points of 
the two lowest BO sheets. We verified that this case 
has no degenerate points anywhere on the JTM, the gap 
between the lowest and next lowest BO sheet being a 

constant -^\qd\ (n)^^^ fld^ across all the JTM. We 
conclude that the P'^' (g d'^^'^ is not a Berry-phase-free 
model. The same applies to all "complete" systems, i.e. 
the pt^^Od''^^^', since {2L 0\L 0,L 0) (corresponding to 
the nondegenerate electronic state) is the largest Clebsch- 
Gordan coefficient among all (2L 0\L M,L — M) = 
N{L)/{L -My.{L + My. (where N{L) is a function of L 
only). 

If we consider now the P'^^^^d^'^^^ models, we see that 
the coefficients \{L 0\L M,L — M)\ are peaked around 
\M\ ^ M Ki 86%L (for large L). For L < 50, say, the 
peak is rather sharp, allowing to concentrate on the 2x2 
block of the interaction matrix: 



-M,-M ^-M,M 











(14) 



Other diagonal and off-diagonal contributions in the 
block vanish, since they could only be related to q^2M^ 
which does not exist, since 2M > I — L. This block, 
therefore, gives a twofold degenerate electronic ground 
state, with second-order departures from the degeneracy, 
as the distortion moves away from the go axis. We con- 
clude that the two lowest BO sheets, in all the 'D'^^'> 1^ d'^^'> 
models, get tangent at (at least) one point, thus making 
equivalent all loops on the JTM. As a consequence, this 
class of models must be considered as Berry-phase-free. 

Finally, we apply the same reasoning to the study of 
a generic V'-^^ (g) d^'^ model. For I < L the \M\ = M, 
at which the relevant coefficients \{l 0\L M,L — M)\ 
is maximum, always satisfies the inequality 2M > I. 
Thus, these are yet more Berry-phase-free systems. On 
the contrary, for I > L, we go from a case without 
(X'(^) (g d(^)) to a case with (I?'^) g) d^^^)) Berry phase, 
thus the result is not trivial. For large enough /, the 
maximum coefficient is attained at some M < I, and 
this leaves nonzero direct off-diagonal coupling elements 
in the 2x2 block ([T^. In such a case, the degenerate 
point is not a tangency, but a conical intersection, thus it 
cannot belong to the JTM. It can be verified numerically 



4 



that the tangency point on the qo axis disappears for all 
I > lc{L), where the "critical value" lc{L) is found be- 
tween L and 2L. For L — 2,3,4, ...10, lc{L) takes values 
4, 4, 6, 8, 8, 10, 12, 12, 14 respectively: thus, for example, a 
V^'^'^ (g) S-^^ model has no Berry phase, while I?*-^-* ® S^^'^ 
can have it. Our demonstration focuses on the go axis: we 
cannot rule out the possibility of other tangency points 
elsewhere on the JTM. Thus the presence of the Berry 
phase is verified only for I = 2L, where the gap between 
the two lowest sheets is a constant on the JTM, while it 
is only likely for lc{L) < I < 2L. 



V. NUMERICAL RESULTS 

We have classified Berry-phase-wise a large class of 
spherical DJT systems: it should be possible to evidence 
the signatures of this property in their vibronic spectrum. 
It was previously shown |^ that the ground state of the 
I?^^-* ® d*-^-* system becomes nondegenerate at strong cou- 
pling. We have applied the same technique of numerical 
diagonalization on a truncated basis to the next Berry- 
phase-free model, I?*-'*-' (g) d^'^'' . In Fig. ||, we show the 
vibronic energy of the lowest L = and L = A vibronic 
states above the BO potential minimum -Ecias^ the low- 
est nondegenerate and degenerate states cross at « 8, 
the former becoming the strong-coupling ground state, 
as the absence of a Berry phase predicts. 

However, an analogous test for the I?*^^' (g) d^^' model 
finds a degenerate L = 4 ground state up to coupling g — 
20. An explanation could be searched in the nature of 
the lowest L = state at weak-coupling: while in V^^^ g) 
niodels it comes as a fragment of the one-vibron 
multiplet, in this case it originates from the two-vibron 
multiplet, due to angular momentum conservation. It 
seems as if the angular moment of the I = 2 vibration 
was insufficient to screen the large electronic moment, 
thus, even in absence of Berry phase, the ground state 
remains degenerate. 



VI. DISCUSSION 

In the Berry-phase-free DJT systems, the tangency of 
the lowest two BO sheets, strictly speaking, invalidates 
the BO treatment, which assumes a large gap between 
the lowest electronic state and the next one. Thus, para- 
doxically, in these systems it is not the Berry phase which 
is related to a breakdown of the BO approximation, but 
its absence. Indeed, we have shown that, even though 
these degeneracies are present only locally on the JTM, 
they affect radically the whole coupled dynamics. 

Our analysis considers spherical DJT models: however, 
it can be extended to molecular point groups, by sub- 
stituting the relevant Clebsch-Gordan coefficients. |1^] 
There are cases, such as the linear T (g) h in icosahedral 



symmetry (or equally coupled T (e + t) in cubic sym- 
metry) which are equivalent to some spherical models, 
but of course this is not always the case, and other sur- 
prising results could be found by looking for tangencies 
in different systems. 

Also, we assume a linear JT coupling scheme (Hamil- 
tonian (|l|)), which is the less realistic, the stronger the JT 
distortion. The perturbative introduction of higher-order 
couplings has usually effects similar to those produced by 
different linear coupling gE ^ gr vc^ cubic symmetry, |lj] 
i.e. of "warping" the JT trough. The continuous JTM be- 
comes a set of isolated minima, connected by low-energy 
paths passing through saddle points. The symmetry of 
the Hamiltonian and of the JTM is reduced to the sym- 
metry group G of the molecule. Yet, the connectedness 
properties are topological properties, therefore they are 
robust against perturbations such as the potential warp- 
ing: even if the symmetry is reduced, the Berry phase 
is still present or absent as determined by the linear 
part. Ham has shown for the E ®e coupling scheme 
that the introduction of quadratic terms in the Hamilto- 
nian does not substantially change the picture as far as 
the Berry phase and the degeneracy of the ground state 
are concerned. In fact, even at strong JT coupling, the 
tunneling among rather deep isolated minima is affected 
by the electronic phase, 1^ and, as a result, the lowest 
tunnel-split state retains the same symmetry and degen- 
eracy as in the purely linear-coupling case. Of course, in 
the extreme limit of very large distortion, higher-order 
terms dominate, and the DJT is replaced by a static dis- 
tortion, where any Berry phase argument becomes irrel- 
evant. 



VII. CONCLUSION 

In summary, we propose evidence for a whole family, 
following the 2?^^-' ® dP'\ of Berry-phase-free dynami- 
cal JT systems: these are the I?^^' ® d^'-* models, with 
I < lc{L), where the critical value lc{L) lies between L 
and 2L. For these models, we also show that the ab- 
sence of the Berry phase does not automatically imply 
a nondegenerate strong-coupling vibronic ground state. 
Moreover, we prove that a Berry phase is present in the 
V'-^^igd^'^^^ models, which, as a consequence, have a 2?^^^ 
[2L + l]-fold degenerate ground state for any coupling. 
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FIG. 1. A sketch of the electronic sphere. The picture in- 
dividuates the two classes of paths mapping onto closed loops 
in the ,JTM: paths of the type 112 involve a sign change (from 
A to A') of the electronic state (a Berry phase). These two 
types of paths can be distinguished in all N > 3-dimensional 
cases. 
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FIG. 2. The energy of the lowest L = and L = 4 
vibronic states of the D^"^^ (g) d''*' system, as a function of 
|£^cias| ~ ^^g^huj. The residual zero-point energy of ^Tiuj is 
subtracted. The energies, in units of the harmonic quantum 
hu!, are obtained by exact diagonalization of the Hamiltonian 
(^) on a truncated Hilbert space including up to 16 boson 
states, enough to reach convergence in this range of couplings. 
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